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Abstract 



Explicit sufficient and necessary conditions for separability of higher 
dimensional quantum systems with rank two density matrices are given. 
A nonseparability inequality is also presented, for the case where one of 
the eigenvectors corresponding to nonzero eigenvalues is a maximally 
entangled state. 
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Quantum entanglement has been investigated for decades because of its importance in 
the foundations of quantum mechanics, particularly in relation with quantum nonsepara- 
bility and the violation of Bell's inequalities [|l|. Recently it has been viewed also as a 
potentially useful resource for communication, information processing and quantum com- 
puting 0, such as for the investigation of quantum teleportation 0, dense coding 0, 
decoherence in quantum computers and the evaluation of quantum cryptographic schemes 
i- 

Due to recent works by Peres |^ and Horodecki et al ^ there exist a simple criterion 
allowing one to judge, whether a given density matrix p, representing a2x2or2x3 
composite system, is separable. Nevertheless, the general problem of finding sufficient and 
necessary conditions for separability in higher dimensions remains open (see e.g. [P, and 
references therein). 

A general condition for separability of a quantum state could in principle be obtained 
from the measure of entanglement. To quantify the degree of entanglement a number of 



entanglement measures, such as the entanglement of formation and distillation [O, O, 13 



negativity [0, Q and relative entropy [|T^, |T^, have been proposed for bipartite states 



[6,8,15-18]. However most proposed measures of entanglement involve extremizations which 



are difficult to handle analytically. For instance, the "entanglement of formation" |Tl 
which is intended to quantify the amount of quantum communication required to create a 
given state, is defined for arbitrary dimension, but so far no explicit analytic formulae for 
entanglement of formation have been found for systems larger than a pair of qubits (spin-^ 
particles). For the case of a pair of qubits, it has been shown that the entanglement of 
formation can be expressed as a monotonically increasing function of the "concurrence" 
variable C. This function ranges from to 1 as C goes from to 1, so that one can take 



the concurrence as a measure of entanglement in its own right ||T^ . 

In this letter, we introduce a generalized concurrence and study sufficient and necessary 
conditions for separability of higher dimensional quantum systems. In particular, we con- 
sider density matrices with rank two. The separability condition for these kind of mixed 
states in arbitrary dimensions is explicitly given. In addition, we present a non separa- 
bility inequality valid in the case where one of the eigenvectors, corresponding to nonzero 



2 



eigenvalues, of a density matrix is a maximally entangled state. 

We first introduce a generalized concurrence C^v for A^-dimensional quantum states. Let 
H be an A^-dimensional complex Hilbert space, with Cj, i = 1,...,N, as an orthonormal 
basis. A general pure state on if ® is of the form, 

N 

1^) = CLijCi^ej, aij G C (1) 

N 

with the normalization ^ = 1 (* denoting complex conjugation). 

Let U denote a unitary transformation on the Hilbert space Ti, such that 

N 

Uci^ YbijCj, bij G C (2) 

N 

and ^2 bijblj = (with Sik the usual Kronecker's symbol). We call a quantity an invariant 
associated with the state |\E') if it is invariant under all local unitary transformations, i.e. all 
maps from H^H to itself of the form U ®U . Let A denote the matrix given by {A)ij = a^j. 
By generalizing the results of analysis on invariants for qubits |]I9|, we can show that the 
following quantities are invariants under local unitary transformations: 

I^=Tr{AA^Y+\ a = 0,1,...,A^- 1; (3) 

(with the adjoint of the matrix A). 

We define the generalized concurrence to be: 



N ^ 2 

^ 2(iV-^T)' ^ I'^ifcCfjm ~ 0'imO'jk\ ■ (4) 



j j',fc,m=l 

For = 2 we have C2 = 2|aiia22 — ai2'^2i|; which is just the definition of concurrence for a 
pure state of two qubits |T^. 

For general iV, we see that when the state |\Ef) is factorizable in the sense that aij = aibj 
for some a^, bj E C, i = 1,...,N, Cn is zero. When |\E^) is maximally entangled, |aij| = 
1/\/N, aij = for z 7^ j, C^v is one. In terms of the Schmidt decomposition, a given |\E') can 
always be written in the form, possibly by changing the orthonormal basis {cj}, i = 1, N: 

N 

1=1 



N 

where ^ Aj = 1, Aj > 0. The invariants are then of the form 

i=l 

N 

4 = 5:A-+\ a^O,...,N-l. 
1=1 

We note that 

AT N 

i=l i^j 
N N 

Therefore = imphes that ^ AjAj = 0. As A, > and ^ Aj = 1, we have that in this 
case only one of the Aj, i = I, ...N is equal to 1 and all other ones are zero. Hence Cn = 
implies that |^) is separable. 

N N 

If Cjv = 1, we have ^ AjA^ = — — — , which is equivalent to the condition ^ A^ = 1/N, 

TV TV 

according to the normalization ^ Aj = 1. The equation ^ A^ = 1/N describes an (iV — 1)- 

1=1 i=i 

N 

dimensional sphere in IR'^ with radius 1/\/N, whereas ^Aj = 1 is a hyperplane in IR^ . 

1=1 

These hypersurfaces have only one contact point at Aj = i = 1,...,N. Therefore 

Cjv = 1 implies that |^) is maximally entangled. We remark that the above properties 
of Cn do not mean that Cat is in general a suitable measure for general A'"-dimensional 
bipartite quantum pure states. It can however be shown that when the matrix AA^ has 
only two different nonzero eigenvalues, the entanglement of formation is a monotonically 
increasing function of Cat, thus Cn can indeed be interpreted as a measure of entanglement 
in this case. 

Let p be a rank two state in H<S)H, with \Ei), \E2) being its two orthonormal eigenvectors 
corresponding to the two nonzero eigenvalues: 

p^p\Ei){Ei\+q\E2){E2\, (5) 

TV 

where q — 1 — p & (0,1). Generally \Ek) = ^ a^^Cj (g) ej, a^j e C, with normalization 

N 

a'lj{a^j)* — 1, k — 1,2. We would like to give an explicit algebraic condition for the 
separability of the above state in terms of af^ . A; = 1, 2. 



With the notations: 



f^ij — ^ij^kl + ^Ij^kl ~ ^il^kj ~ ^il^kj^ 



we have the following conclusion: 

[Theorem 1]. If all a'^j are real, p is separable if and only if one of the following quantities 
(Ai or A2) is zero : 

ijkl ijklmn 

A2 = El7S + (i-p->gr + El/^r; 

ijkl ijkl 

or equivalently one of the following two sets of relations ((6) or (7)) hold: 

7g = (l-p->g' PiKn-o^lPmn, yi,j,k,l,m,n. (6) 
7g = -(l-p->g, /3g = 0, yt,j,k,l. (7) 

N 

[Proof]. From the discussions above we have that a state |^) = E o^ij^i^^j is separable 
if and only if the generalized concurrence Cjv is zero, i.e., 

aijdki ^ aaakj, \/i,j,k,l. (8) 

Therefore a vector of the form + A|£^2), A e C is separable if and only if A is a common 
root of the following equation set Eqij-kf. 

+ /^S"^ + 7g = 0, yt,j,k,L (9) 

I 

We first prove the necessity part of the theorem. Suppose that p — '^p[\Ut){Ut\, with 

t=i 

I some positive integer and 0<p^<l,X^pj = l, gives a decomposition of p in terms of 
separable normalized vectors Ut (where \Ut){Ut\ means the projection onto the vector Ut). 
As each of these vectors must be in the range of p (since for each fixed t, \Ut){Ut\ < p[~^p), 
we can write them as linear combinations of the two eigenvectors and \E2) which span 
the range of p. So let us assume that Ut — c\\Ei) + cl\E2) (for some c\,cl e C). We deal 
with the problem in two cases: 

Case 1. I £'2) is not separable 

In this case, since c\ can not be 0, Ei + \tE2 will be a separable vector, where \t = c\/c\. 
Clearly, not all s can be equal, otherwise all the Ut s would be constant multiples of a 
fixed vector, and p would be of rank 1. Hence there must be at least two distinct choices of 
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A(. On the other hand, as l-E'2) is not separable, its generahzed concurrence Cat is not zero. 
Hence there is some io,jo,ko,lo such that 7^ 0, i.e. the relation Eqi^j^-k^i^ is indeed a 
quadratic equation. It must have exactly two roots, say A^^\ A*^^\ and so these two values are 
the only possible choices for the At's. But in order that there is not only one possible choice, 
the above two roots must be different. And all the relations Eqij.^i have these different two 
roots. Consider for any j, k, /, 

Situation 1: If afj 7^ 0, the corresponding relation is then not an identity. All the 
quadratic equations in the equation set Eqij-^i have the same two distinct roots. In other 
words, Eqij-ki and Eqi^j^.^Qi^ have the same roots. From the standard theory of quadratic 
equations, we have 

i-^tj 1030 r^*ojo «j ' V-^"/ 
fcZ fcoio ^ rykolo^kl Ql N 

Situation 2 : If afj = 0, then the equations Eqij.^i become identities, i.e. P^j and j^j must 
be too, because otherwise at least one of the relations Eqij.^i would be a linear equation. 



and so could not be satisfied by two distinct roots. Thus in this case (JWj) and ([IT]) also hold 
(as both sides in both equations vanish). 

Thus we have obtained ([l0|) and (|lT]) for arbitrary ijkl and for all ioJo^o^o such that 
af°j° 7^ 0. As (|10D and (|lT]) are trivially valid whenever ijkl and iojokolo are such that 
a^j = af°j° = 0, we have indeed established ([l0| ) and (pT]) for all quadruplets ijkl, ^oJo^o^o- 

We now denote the above two distinct roots, which are common to all of the equations 
Eqij-ki, by /ii,/i2, with the convention that in case the roots are real, ^2 will denote the 
larger one. Each vector Ut is either of the form Ei = -^1+^^-^^ or of the form E2 = ^i±MiM3=. 
Therefore we can write p as, p = p'\Ei){E2\ + (1 — p')\E2){E2\, with < p' < 1. Comparing 
with the coefficients of \Ek){Ei\, k, I = 1,2 in the expression (|^), we get that the above 
decomposition of p is equivalent to the following two relations: 

+ 1,1.. 12 =P^ 1 , I.. 12 + 1 , I.. 12 = 0- (12) 



l + |pi|2 I + IP2I' 1 + l + W^ 

Solving (|T^) for p and p' we get 



p=(l-piP2-) , p' = — , (13) 

z z — P1P2-2 



where ^ = /i2 — /ii- 

As afjS are real numbers, fii, fi2 are roots of a quadratic equation with real coefficients. 
For a quadratic equation ax'^ + bx + c = with a,b,c & M, a ^ 0, and roots a, (3 with 
a 7^ /?, we have that is either +1 or —1, depending on whether 6^ — 4ac is positive or 
negative respectively. Hence in our case we have that fiifi2 = +(1 ~ P~^) or —(1 — p~^). 
On the other hand, since fiifi2 is real, the solution for p' in ( pTSD implies that ^^^3^ is real, 
which is possible if and only if either the roots are both real or the roots are both purely 
imaginary. 

For the the case where the roots are both purely imaginary, we have, yUi/i2 = —(1 —p^^)] 
and a direct simplification from (|T3p gives p' = ^. The condition for having purely imaginary 
roots of quadratic equations gives that /?g = 0, ^ijkl. 

In the case where the roots are real, /i2 being the larger root by assumption, we have 
fJ'ifJ'2 = (1 — P^^)- By replacing by in the expression of p' in (0), we get that the 
condition that p' takes values between and 1 is equivalent to > 0, /xi < 0, which is 
trivially valid since fi2 > /ii and /ii/i2 < in this case. 

It remains to observe case that yUi/X2 is nothing but the ratio which is either 1 —p^^ 
or — (1— p^^). Taking into account relation (pT]), one concludes that either 'j^j = (1— p^^)a;fj 
or jfj = —(1 — p~^)afj for all i,j, k, I. It then follows that in the second case jS^j = for all 

k, I. This completes the proof for the necessity part in the case 1. 

Case 2. I-E2) is separable 

In this case from (||) we have a^j = 0, Vi, j, fc, /. Since not all of the f/t's can be multiples 
of l-E'2), we must have at least one choice of A such that Ei + A-E2 is separable. This A 
must be a common root to all equations Eqij-ki as before. All these equations are now 
linear ones. Excluding the trivial case that all of them are identically (which means that 
\Ei) is also separable), we see that there is only one possible choice of A (which is the 
solution to a nontrivial linear equation). Then p can be expressed as p = p"\E2){E2 \ + (1 — 
p") l-^i+'^-^-O^C^i+'^'^al ^ ^i^g^i ig p" = X, which is clearly a contradiction. Thus in case where 
\E2) is separable, we must have \Ei) to be separable too in order for p to be separable. It 
is clear that in this case the conditions of the theorem hold. 
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Now we prove the sufficiency part for the theorem. If either (Q) or (|^) holds, from the 
previous analysis it is clear that the equations Eqij-^i have common roots. If |i?2) is not 
separable, then not all of these equations are identities. And there are at most two common 
roots. Furthermore, if holds, the product of the two roots must be (1 —p~^) < 0, so that 
the two roots be real and unequal, with opposite signs. Denoting them by /ii, ^2 as before, 
we can easily get from the arguments given earlier an explicit decomposition of p in terms 
of separable pure states. Similarly, if (|^) holds, the two roots must be purely imaginary 
(since /S^j = 0) and the rest also follows easily from our earlier analysis. Finally, if \E2) is 
separable, either of the conditions (||) or forces •y^j to be for all ijkl. Hence \Ei) is 
separable too and p is trivially separable. This completes the proof. 

The conclusions in Theorem 1 can easily be generalized to the complex case, af^ G C 

[Theorem 2]. p is separable if and only if there is 6 & M such that 

7g = e^'(l-p"'X (14) 
= , yt,J,k,l,m,n; (15) 

and 

^^^ii±^G[0,l], (16) 
z - P1P2Z 

where z = e*^z, z = p2—pi 7^ 0, pi and p2 are the roots of the equation a'^jX^+P^jX+'j^j = 0, 
for some k, I such that a^j 7^ 0. 

The proof for the necessity part of the theorem is similar to the proof of the corresponding 
part in Theorem 1. One only needs to note that since z/z is of modulus 1, and of the form 
some 9. From the relations (P^Sl), we obtain then that ^ is a real number 

between and 1. 

The proof of the sufficiency part for the theorem essentially follows by observing that 
the arguments given in the necessity part can be reversed. From the previous analysis, it 



follows from ([M]) and ([T5|) that the equations Eqij-^i have common roots. In case all 



a. 



kl 



are 0, by the condition 'yfj = e*^(l — p ^)afj, we see that all 'yfj are too. Hence both 

,.khr 



\Ei) and |i?2) are separable, and so is p. In the other situations where some of the afj's 
are nonzero, the corresponding equations Eqij-^i have exactly two roots which are different 
by the condition (|1^). Moreover, (p!6| ) ensures that the explicit solution for the choice of 



p', which gives the decomposition of p in terms of separable states, is indeed a vahd choice 
between and 1. This completes the proof of theorem 2. 

TV 

[Corollary]. Let \E2) be the maximally entangled vector given by \E2) = ® 
For any vector \Ei) which is orthogonal to 1-^2); P = p\Ei){Ei\ + (1 — p)\E2){E2\ is not 
separable for < p < |. 

[Proof]. From the necessary condition for separability, 7^- = e*^(l —p~^)afj, we get that 



= i^C(2), where 



N 



TV 

E 



^ ^ are 



^ 2(7V3T) I7ijl diiLi o(2) = 2^ Fij- 

\ i,j,k,m=l \ i,j,k,m=l 

the generalized concurrences associated with the states \Ei) and I-E2) respectively. As a 
pure state is separable if and only if the corresponding generalized concurrence is zero, a 
necessary condition for separability is that C(i) and C(2) should be inversely proportional to 
the contribution of the corresponding pure state to p, i.e. the eigenvalue corresponding to 
that eigenstate. From = ^— ^ < 1, we have p > h. 

We have studied the sufficient and necessary conditions for separability of higher di- 
mensional quantum systems of density matrices with rank two, in terms of the generalized 
concurrence. Our approach concerns the eigenstates of a given density matrix p (thus results 
as in our corollary are straightforward to deduce). It gives a direct minimal decomposition 
of separable states in terms of pure product states. The method can also be used to check 
whether it is possible to construct a separable mixed state from two given pure states. 

Another elegant approach, based on the positive partial transpose (PPT) condition, deal- 
ing with the sufficient and necessary criteria for separability of higher dimensional quantum 
systems was presented in ||2^ 0. For a given density matrix p, one can judge its separabilty 
by calculating its supporting dimensions (p is supported on M x if this is the smallest 
product Hilbert space on which p can act), checking the positiveness of partial transposed 
p and comparing the rank of p with dimensions M, A^. 
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